Availability and readiness of healthcare facilities and their effects on long-acting modern contraceptive use in Bangladesh: analysis of linked data

Aim Increasing access to long-acting modern contraceptives (LMAC) is one of the key factors in preventing unintended pregnancy and protecting women’s health rights. However, the availability and accessibility of health facilities and their impacts on LAMC utilisation (implant, intrauterine devices, sterilisation) in low- and middle-income countries is an understudied topic. This study aimed to examine the association between the availability and readiness of health facilities and the use of LAMC in Bangladesh. Methods In this survey study, we linked the 2017/18 Bangladesh Demographic and Health Survey data with the 2017 Bangladesh Health Facility Survey data using the administrative-boundary linkage method. Mixed-effect multilevel logistic regressions were conducted. The sample comprised 10,938 married women of 15–49 years age range who were fertile but did not desire a child within 2 years of the date of survey. The outcome variable was the current use of LAMC (yes, no), and the explanatory variables were health facility-, individual-, household- and community-level factors. Results Nearly 34% of participants used LAMCs with significant variations across areas in Bangladesh. The average scores of the health facility management and health facility infrastructure were 0.79 and 0.83, respectively. Of the facilities where LAMCs were available, 69% of them were functional and ready to provide LAMCs to the respondents. The increase in scores for the management (adjusted odds ratio (aOR), 1.59; 95% CI, 1.21–2.42) and infrastructure (aOR, 1.44; 95% CI, 1.01–1.69) of health facilities was positively associated with the overall uptake of LAMC. For per unit increase in the availability and readiness scores to provide LAMC at the nearest health facilities, the aORs for women to report using LAMC were 2.16 (95% CI, 1.18–3.21) and 1.74 (95% CI, 1.15–3.20), respectively. A nearly 27% decline in the likelihood of LAMC uptake was observed for every kilometre increase in the average regional-level distance between women’s homes and the nearest health facilities. Conclusion The proximity of health facilities and their improved management, infrastructure, and readiness to provide LAMCs to women significantly increase their uptake. Policies and programs should prioritise improving health facility readiness to increase LAMC uptake.


Introduction
Unintended pregnancy represents around 44% of all pregnancies in low-and middle-income countries (LMICs) [1]. Almost one-third of them occur due to contraceptive failures, primarily because of using short-term modern contraceptives (e.g., pills and condoms) and traditional contraceptives (e.g., withdrawal, rhythm) [2]. In LMICs, unintended pregnancies are responsible for around 55 million unplanned births and 25 million miscarriages in a year [1]. These constitute around 80% burden of pregnancy complications (including severe bleeding and infections) [3], leading to over 75% of the 118,000 maternal deaths each year in LMICs [4]. The underlying reasons are adverse health and behavioural issues following unintended conception, during pregnancy and the postpartum period, including depression, anxiety and lower use of intra-partum, delivery and post-partum care [5][6][7]. Furthermore, such high occurrences of unintended pregnancies in LMICs are attributed to 138 million abortions annually [1]. Abortion is prohibited in many LMICs, even though having an abortion can save a woman's life. This prohibition and scarce safe abortion facilities result in unsafe abortions, increasing maternal morbidity and mortality. Also, unwanted pregnancy forces many women to carry their fetuses to term and postpartum and experience the socioeconomic strain of the subsequent births and health risks. Access to a full range of contraceptives and family planning services means that women are more likely to find something that works in their circumstances and increase the likelihood of consistent and correct use. Indeed, contraception use can avert about 60 million healthy lives lost [8] and around 10% of child deaths in LMICs per year [9]. The injectable, female strelization, male strelization, intrauterine devices (IUD) and hormonal implants are the long-acting modern contraception (LAMC) that can significantly reduce these burdens [10]. Accordingly, universal access to sexual and reproductive healthcare services, including increasing use of LAMC, has been targeted in the Sustainable Development Goals (SDGs, Targets 3.7) to be achieved by 2030 [11]. This is also considered an essential pathway to achieving another two SDG targets: reducing maternal (Target 3.1) and child mortality (Target 3.2).
In Bangladesh, an LMIC, around half of the pregnancies are unintended at conception, and over one-third are attributed to contraceptive methods failure [2,12]. Consequently, the government has focused on increasing contraception coverage, particularly LAMC use [12]. Contraceptives are being provided from all existing healthcare facilities and community clinics, strengthening the traditional approach of providing family planning services and contraception by family planning workers [13]. Targeted programs are also implemented to increase health facilities' capabilities and readiness to offer a full range of contraceptive methods [12]. The target is to enable couples to access contraceptive methods that meet their needs. A similar approach has been taken in many other LMICs. However, to what extent this endeavour contributes to an increase in LAMC use is an ongoing debate in Bangladesh, as the prevalence of modern contraceptive use has been minimal over the years, along with evidence of rising pregnancies due to contraceptive failures [13]. As a result, policymakers are in a dilemma about the contribution of healthcare facilities to LAMC use and cannot decide whether they should prioritise the health facility over the traditional approach of family planning services at the community level by family planning workers [14]. This dilemma, to some extent, is due to insufficient evidence on how contraception availability in health facilities affect contraceptive use -a lack that also suppresses the true effects of population-level factors [15][16][17][18]. Thus far, limited research has been conducted mostly in African countries to examine the effect of health facility-level factors on the overall use of contraception [19][20][21][22]. However, those studies did not present a stratified analysis for rural and urban areas, despite substantial variations in the distribution of health facilities and contraception-using patterns in these two settings. The norms of contraceptive use and associated factors for Bangladesh are significantly different from African countries, along with a substantial rural-urban variation [2]. However, previous studies on this topic in Bangladesh considered only the population-level factors [2,23,24]. To address these gaps, we conducted this study to examine the effects of availability and accessibility of health facility-level factors on LAMC use in Bangladesh.

Study overview
We analysed data extracted from 2017/18 Bangladesh Demographic and Health Survey (BDHS) and 2017 Bangladesh Health Facility Survey (BHFS). We established linkage between these two survey datasets using the administrative boundary linkage method [25,26]. The Demographic and Health Survey Programme of the USA designed these surveys, and the National Institute Keywords: Modern contraceptives use, contraception, linked data, health facility data, multilevel regression, Bangladesh of Population Research and Training conducted the surveys under the supervision of the Ministry of Health and Family Welfare of Bangladesh. Detailed sampling procedures of these surveys can be found in the respective survey reports [27,28]. Briefly, the 2017/18 BDHS collected nationally representative sample selected by using multistage stratified random sampling methods. At the first stage of sampling, 675 clusters were selected from a list of 293,579 clusters that was created as part of the 2011 National Population Census. A fixed number of 30 households were selected at the second stage from each cluster through probability proportional to the unit size. Of them, data were collected from 19,457 households. There were 20,376 eligible women in the selected households of which 20,127 women were included in the survey with a response rate 98.8%. The 2017 BHFS used a list of 19,811 registered health facilities as a sampling frame generated by the Ministry of Health and Family Welfare. A total of 1600 health facilities were selected from this list, and 1524 health facilities were finally included in the survey.

Sample
We analysed data from 10,384 women selected from 672 clusters included in the BDHS 2017/18. The criteria used for inclusion in this study were women of reproductive age , who were fertile but not pregnant or experiencing lactational amenorrhea at the time of the survey and did not desire a child for 2 years prior to the date of the data collection.

Outcome variable
The outcome variable was LAMC use. The relevant data were collected by asking women, "Are you currently doing something or using any method to delay or avoid getting pregnant?" Women who responded "Yes" to this question were then asked, "Which method are you using?" and could select their method from the following list: pill, injectable, implant, intrauterine device (IUD), condoms, female sterilisation, male sterilisation, periodic abstinence and withdrawal. A free text option was also provided to report the name of contraception if not included in the list. If a woman reported multiple methods, the method they used most frequently was selected. Responses were coded as "1" for LAMC (i.e., the injectable, implant, IUD, female sterilisation and male sterilisation) and "0" for contraceptive non-use (responded "no" to the first question) or for use of traditional methods (i.e., periodic abstinence or withdrawal). A further classification was made as use or non-use of long-acting reversible methods (injectable, IUD and implants) and parament methods (male sterilisation, female sterilisation). Women who reported condoms or pills as their current method of contraception were excluded, as these are available at local shops and pharmacies in Bangladesh as well as supplied by the family planning staff working at the field level.

Explanatory variables
The explanatory variables in this study were selected in three stages: at the first stage, we conducted a comprehensive literature search in five databases (PubMed, CIN-HAL, Web of Science, Embase and Google Scholar) and identified all relevant studies conducted in LMICs, particularly in Bangladesh [15,16,18,23,24,[29][30][31][32][33][34][35][36]. The variables considered in these selected studies were then summarised, and their availability was then checked with the dataset we analysed in the second stage. The summarized variables that were available in the survey were then considered with the outcome variable of this study and forward regression models were run. Less significant variables (p > 0.20) were removed from the list. Furthermore, we checked multicollinearity before entering these variables into the model. If evidence of multicollinearity was found (VIF > 10), the relevant variable was also deleted. Finally, the selected variables were classified into four broad groups in line with the multilevel analytic approach: health facility-, individuals-, household-, and community-level factors.
The health facility-level variables were general service readiness (facility's management and infrastructure), LAMC services availability and readiness scores. For facilities that provide contraceptions, indices of severity of the availability of LAMC at the nearest health care facility and readiness were created following the WHO guidelines [37]. Seven variables were used to create contraception availability scores: combined oral contraceptive pills, progestin-only contraceptive pills, progestin-only injectable contraceptives, intrauterine device, implant, male sterilisation and female sterilisation. These indicator variables were classified in line with this study's aim to generate LAMC availability scores: 1 point for each form of LAMC availability in the health facility and 0 for non-availability. Similarly, the contraceptive service readiness score was generated using seven dichotomous variables of healthcare personnel availability and training on providing LAMC. These scores were generated using the principal component analysis. The average distance on-road communication of the nearest LAMC-providing facility was also calculated and included as a health facility-level variable. Cluster's LAMC-providing nearest health facility was identified first. Using Bangladesh road communication data, the average distance from the cluster to the nearest health facility was calculated separately for each of the eight administrative divisions. We used a regional-level average instead of the actual distance between the clusters and closest health facilities because the BHFS survey included a sample for all health facilities except the district health facility and maternal and child welfare centre. Thus, a cluster's nearest health facility might not have been selected and included in the survey; hence, the actual distance would be problematic. The details of these computation procedures can be found elsewhere [19].
The individual-level variables are participants' age (≤19 years, 20-34 years, ≥35 years), education (no formal education, primary, secondary, higher), and employment status (yes, no). The household-level variables are partners' education (no formal education, primary, secondary, higher), partners' occupation (agricultural worker, physical worker, service, business, others), household wealth quintile (poorest, poorer, middle, richer, richest) and the number of children ever given birth (no child, 1-2 children, > 2 children). BDHS generated wealth quintile variables using principal component analysis of the relevant data related to household assets. The scores generated were divided into five equal groups with a cut-off value of 0.20 and consecutively defined as poorest (0.20) to richest (1.00). Family types were also included as a household-level variable and classified as Nuclear family (if the family had < 5 members) and joint family (if the family had ≥5 members). The community-level factors were the place of residence (rural and urban) and administrative divisions.

Statistical analysis
Descriptive statistics were used to summarise the characteristics of the respondents. Multilevel logistic regression was used to assess the associations of LAMC use with health facility-, individual-and household-, and community-level factors. The reason for using the multilevel regression was the hierarchical structure of the BDHS data, in which individuals are nested within a household and households are nested within a cluster. This created multiple dependencies in the data for which multilevel regression was deemed appropriate [38]. We followed the progressive model-building technique and developed four different models. Model 1 was the null model where no covariate was adjusted. Model 2 was the health facility-level model where LAMC use was considered with the health facility-level variables. Model 3 was the extension of Model 2 with the inclusion of individualand household-level factors. All health facility-, individual-, households-, and community-level variables were adjusted in Model 4, which was the final model. Results were reported as adjusted Odds Ratios (aOR) with 95% Confidence Intervals (95% CI). The Intra-Class Correlation (ICC), Variance Inflation Factor (VIF), Akaike's Information Criteria (AIC) and Bayesian Information Criteria (BIC) for each model were recorded and compared to select the best model. The ICC value was calculated by dividing the between-clusters-variance of LAMC use (random intercept variance) by the total variance of LAMC use (sum of between-clusters-variance and within-cluster (residual) variance of LAMC use). Statistical package R was used for analyses. All methods were performed in accordance with the relevant guidelines and regulations. Table 1 shows the background characteristics of the respondents. A majority of the respondents in this study were in 20-34 years age group (75%) at the time of the survey. Almost a quarter of women received no formal education, a third completed only primary and another third completed secondary education. Approximately 4% of women had no children and over 60% had 4 years or more intervals in their two most recent pregnancies. Around 74% of women resided in rural areas.

Distribution of long-acting modern contraceptives in Bangladesh
The distribution of long-acting modern contraceptives from health facilities in Bangladesh is presented in Table 1. We analysed data of 10,384 women, 3483 (34%) of them used LAMC. Injectable was the dominant LAMC, used by 55% of the overall LAMC users. Female sterilisation (25%), male sterilisation (6%), and implants (11%) were the other three major forms used by over 42% of LAMC users.

Distribution of long-acting modern contraception use across women's socio-economic characteristics
The overall and method-wise distribution of LAMCs are also presented in Table 1. Of the women who used LAMC, a majority of them were 20-34 years old (79%) and received either primary (40%) or secondary education (32%). Women engaged in formal income-generating activities (60%) were more likely than others to report using LAMC. Around 58% of the LAMC users had more than two children during the survey. The prevalence of LAMC users was higher among women in the Dhaka division (23%) than Rangpur and Chattogram divisions (16%). A similar pattern was recorded for the methodwise distribution of LAMCs. Male and female sterilisations were more common among people with two or more children, compared to IUDs, implants and injectables. Similarly, sterilisations were higher among those who reported having > 4 years of intervals between the two most recent live births (Table 1).

Distribution of health facilities in Bangladesh
The distribution of health facilities across divisions is presented in Table 2. Of the 1524 health facilities, 1357 (89%) provided LAMC, and 93% (1262/1357) of these facilities were government hospitals/clinics located at the Upazila (i.e., sub-district) or in communities. Relatively high proportions of these facilities were located in the Dhaka (n = 245) and Chattogram (n = 234) divisions. Sylhet division had a lower number of health facilities with only 87 Upazila/community level governmental hospitals. NGO clinics (n = 19) and private hospitals (n = 9) were found highly clustered in the Dhaka division. The average distance between the health facilities and the BDHS clusters was around 6 km, higher in the Sylhet division (8 km) and lower in the Dhaka division (5 km). The mean score for the health facility management system was 0.79, higher in the Dhaka division (0.93) than in Khulna (0.87) and Chattogram (0.82) division. Similarly, higher availability of different LAMC was reported for Dhaka (mean score, 0.91) and Chattogram division (mean score, 0.74) than the corresponding average score of 0.69. Of the facilities where LAMCs were found available, 69% of them were functional and ready to provide LAMCs to the respondents. A higher score was reported for the healthcare facility of the Rangpur division (0.88) than for the Dhaka (0.85) and Chattogram (0.86) divisions.

Model selection
The associations of LAMC use with health facility-, individual-, household-, and community-level variables were assessed through multilevel logistic regression models. Of the four models ran separately, the AIC, BIC and ICC values were compared to identify the best models. The preferred model was the one that had the smallest AIC, BIC and ICC (Table 3). According to these markers, Model 4 (including health facility-, individual-, household-, and community-variables) fitted the data better. The ICC value for the null model (Model 1) suggested around 17% differences in LAMC use across clusters included in the BDHS. However, this value was reduced to only 3% once health facility-, individual-, household-, and community-level factors were included in the final model. Around 11% of such reduction was reported once health facility-level factors were included in the null model.

Factors associated with long-acting modern contraception use
In the final model, we found that health facility-level factors were strong predictors of LAMC use following the adjustment of health facility-, individual-, household-, and community-level factors ( Women engaged in formal employments and with increased years of education were more likely to report using LAMC than women who received no formal education and were not employed in a formal job. Women who received primary, secondary, or further education and whose partners received secondary or higher education were significantly more likely than   Table 3 Multilevel logistic regression models assessing the relationships between the use of long-acting modern contraception and health facility-, individual-, household-, and community-level factors (N = 10,384)

Effects of health facility environment on long-acting modern contraception use
We calculated the health facility environment by considering all health facilities providing LAMC in 6.3 km (division level average distance on road-communication of health facility from the BDHS cluster) buffer distance. The number of health facilities and their characteristics were generated for national and urban and rural locations ( Table 4). The likelihood of LAMC use was twofold among women living in areas where there was a health facility within 6.3 km (aOR 2.01; 95% CI 1.67-2.37). The likelihood increased further (aOR 3.39; 95% CI 1.84-4.90) if the 6.3 km buffer distance had more than one health facility compared to having no health facility. These likelihoods increased even further for rural areas. The general health services readiness, degree of availability of LAMCs at the nearest health care facility, and readiness of the nearest health care facility (where LAMCs are available) to provide LAMCs were also found significantly associated with LAMC use (Table 4).
We also conducted a sub-group analysis to explore the association between health facility level factors with permanent modern methods and long-acting reversible methods used (results now shown in the table). The observed associations were found consistent with the observations reported in Tables 3 and 4. However, the associations were stronger for permanent than reversible ref Reference group, ** p < 0.01, * p < 0.05

Discussion
This study provides new insights for LMICs and Bangladesh on the role of availability and accessibility of health facilities in determining LAMC use. Only 34% of the total respondents analysed used LAMC. LAMC methods are available in 89% of the health facilities in Bangladesh; however, 65% of them are found functionally ready to provide these services. The findings suggest the average distance of a health facility that provides LAMC from the respondents' homes was 6.36 km, highest in Sylhet division (8 km) and lowest in Dhaka division (5 km). The LAMC uptake was positively associated with LAMC availability at the nearest health facility and improved management and infrastructure of those facilities. These relationships were stronger in rural than urban areas. The number of health facilities providing LAMC within 6.3 km was revealed as the strongest predictor of LAMC use. Women from lower socio-economic backgrounds were more likely than women of other wealth groups to use LAMC. Overall, our findings suggest that LAMC uptake would increase among women if there are available LAMC-providing facilities with good infrastructure in proximity to women's homes. Although the SDGs target to ensure universal coverage of sexual and reproductive healthcare services by 2030 [39], the Bangladesh government has set the target to achieve 75% coverage of modern contraceptive use by 2025 [40,41]. Achieving this target is deemed a real challenge as only 54% of the eligible couples in Bangladesh use modern contraception [27], and this percentage is even lower for LAMCs (34%), as found in this study. This estimated prevalence of LAMCs in Bangladesh is lower than the global (47.6%) and LMICs prevalences (47.0%) but higher than in the African region (19%) [42]. The prevalence of LAMCs was also reported higher in two neighbouring countries, India (41%) and Myanmar (38%) [15].
In Bangladesh, the programmes for improving health facility management, infrastructure and accessibility to LAMCs were first taken into account as part of the 7th Five-Year Plan (2016-2020) [43] -and are currently functional as part of the 8th Five-Year Plan (2021-2026) [44]. However, no notable progress has been observed to date. Providing contraception services, including LAMC, via family planning workers (who visit couples' homes fortnightly) is still the fundamental way to provide family planning services and contraception at the community level [27]. An advantage of this approach is that couples who are not interested in accessing services from health facilities can still access contraception services. However, this approach has two major disadvantages: i) family planning workers usually lack specialised medical education and ii) family planning workers are restricted to provide non-surgical contraception methods, mainly pills and condoms [2]. Consequently, pills (25%) and condoms (7%) are the main contraceptives in Bangladesh and are used by 60% of the overall contraception users [27]. This pattern differs from many other LMICs where condoms and female sterilisation are the dominant contraception methods [13] These indicate the dominant role of family Table 4 Multilevel logistics regression assessing the relationships between of LAMC use and attributes of health facilities located within 6.36 km distance from the BDHS clusters a Models adjusted with women's age, education, employment status, number of children ever given birth, Intervals between the two most recent live births, husbands' education, occupation, family type, household wealth quintile, place of residence and region ** p < 0.01, * p < 0.05

Health facility characteristics
Long-acting modern contraception use, aOR (95% CI) planning workers in providing contraception services over health facilities [2,28]. Furthermore, the failure rates of pills (5-15%) and condoms (7-19%) are very high in their typical use [46]. Therefore, the ongoing provision of contraception services via family planning workers needs to be improved to reduce unintended pregnancies and related adverse consequences [2] as well as to achieve SDGs targets to reduce maternal and child mortality.
Although several programs to increase LAMC use have been introduced in Bangladesh, including payment for sterilisation or vasectomy, progress is unlikely to be significant unless the provider-level challenges are addressed [47,48]. The first challenge comes from the lack of coordination between family planning workers and health facilities. In Bangladesh, family planning workers work under the Directorate General of Family Planning Services, and health facilities function under the Directorate General of Health Services [49]. One of the job responsibilities of family planning workers is promoting awareness about LAMC use among women in rural areas. They educate eligible couples on family planning methods and contraception use. They do not provide LAMC to interested people but refer them to health facilities [50]. As a result, the coordination between family planning workers and health facilities is minimal [2]. Although family planning workers have direct contact with eligible couples, they cannot monitor whether they accessed health facilities for LAMC or not. Such a lack of coordination results in unequal coverage and a significant pathway of miss-out from the services [13]. Effective linkages between family planning workers and health facilities are crucial. However, because of family planning workers' direct access to the grassroots level, LAMC use among disadvantaged groups and in rural Bangladesh is relatively high, which is unique to the results of studies conducted in LMICs, including Bangladesh [16,17,32,33]. However, LAMC use is still relatively low, even among the women of the poorest and poorer wealth quintiles.
The second challenge comes from the health facility. Because of the government's efforts, LAMC is available in most healthcare facilities in Bangladesh. As per the results of this study, it is available in 89% of healthcare facilities. However, 69% of them are functionally ready to provide LAMC with high divisional variations. Over 85% of health facilities in Rangpur, Dhaka, and Chattogram divisions had functional readiness to provide LAMC. Nearly half of the healthcare facilities in the remaining five divisions were not functionally ready to provide LAMC, although services were available there. This pattern has created several challenges at the field level, including the crowdedness in the nearest LAMC-providing healthcare facility and increased distance between women's homes and that facility.
Availability of LAMC-providing health facilities in short distances has been identified as a critical predictor of LAMC uptake in this study, a consistent finding in other LMICs such as Turkey and Zambia [34,36,51]. This finding also indicates that some people may not access LAMC due to the relatively high distance of the LAMC-providing facilities despite having an intention. There may have a compound effect of this disadvantage among women with low educational achievement and relatively low knowledge of LAMC. Social barriers and taboos, a low agency of decision making and financial inability can prevent less educated women from visiting health facilities to access LAMC if located far away [2]. However, this understanding may not be true for younger women, who are comparatively educated and enjoy more autonomy than older women. Moreover, limited and distant LAMC-providing healthcare facilities are often overcrowded. Furthermore, some healthcare facilities are not women-friendly; for instance, they do not have private space to discuss LAMC use [2]. These challenges create inequality in LAMC use, wherein LAMC-providing healthcare facilities (> 93% are government owned) tend to serve relatively poor people. On the other hand, educated and wealthy people choose private clinics to access LAMC (which are limited in Bangladesh) or use short-acting methods, including condoms or pills from private sources. These observations highlight the requirement for more healthcare facilities and infrastructural development of the existing healthcare facilities that can provide LAMC at the community level. Engaging 13,221 community clinics that are currently functioning could accelerate the initiative.
Our findings also suggest improvement in the management of health facilities has significant contributions to the increased use of LAMC. In Bangladesh, there is a significant shortage of trained healthcare professionals assigned to LAMC and other contraception services [52]. Moreover, healthcare providers are expected to provide counselling in conjunction with LAMC [44]. However, many practitioners, particularly those who completed a Bachelor of Medicine or Bachelor of Surgery degree, may not prioritise counselling over other competing tasks. Moreover, in most healthcare facilities, a substantial proportion of healthcare providers are male, and women usually do not feel comfortable accessing contraception services through a male provider [52]. These challenges are even deep-rooted in rural health facilities along with poor management, including inadequate medical supplies, inability to use modern technologies in providing health services and absence of assigned healthcare personnel [53]. Consequently, women in rural areas are less likely to use LAMC.
This study has several limitations. Data analysed in this study were extracted from cross-sectional surveys. Therefore, the findings are correlational only, not casual. To secure the privacy of the respondents, BDHS displaced clusters' locations up 0-5 km for rural areas and 0-2 km for urban areas. Therefore, the calculated average distance of the nearest health facility could be slightly different from the actual distance. However, the BDHS ensured that the new perturbed locations fell within the designated administrative boundaries. Therefore, errors from displacement are likely to be random and minimum. A previous study found that the effect of this variation is insignificant [54]. Moreover, the contraceptive methods that were included to generate LAMC variables widely vary from one another, and their seeking processes are also different. For instance, the female sterilisation methodd need to be performed in -a well-developed in-patient settings by expert healthcare personnel immediately following childbirth. On the contrary, injectable methods are comparatively easy, do not need a high-end in-patient settings and highly skilled healthcare personnel, and can be applied at any time. These disparities may limit the interpretability of this study's findings. However, we have conducted a sub-group analysis to determine the association of health facility level factors on permanent modern and long-acting reversible modern methods. The association reported were consistent with our main study findings. This confirms that the findings of this study are still valid. Regardless of these limitations, to our knowledge, this is the first study in the context of LMICs that examined the influence of health facilitylevel factors adjusted with the individual-, household-, and community-level factors. Clusters of LAMC use and non-use were identified, and the factors associated with LAMC were determined using multilevel regressions with adjustments for a wide range of factors. Therefore, the finding of this study is likely to be precise and can be generalisable in countries with similar features to Bangladesh and may help develop evidencebased policies.

Conclusion
During 2017-2018, around 33% of women used LAMC in Bangladesh. Health facility-level factors, including health facility infrastructure, management, and readiness to provide LAMC were strong predictors of LAMC use in Bangladesh and the effects of the predictors were substantial in rural than urban areas. Health facilities providing LAMC located at a closer distance were positively associated with increased uptake of LAMC. Health facilities should be strengthened to provide LAMC and more health facilities of such capacities are needed. The current provision of family planning services through the field-level workers also needs to be strengthened with strong coordination and referral linkages with the health facilities.